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"X-Ray Edge" Singularities in Nanotubes and Quantum Wires with Multiple 

Subbands 
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Band theory predicts a (e — A)~^" van Hove singularity in the tunneUng density of states at 
the minimum energy (A) of an unoccupied subband in a one-dimensional quantum wire. With 
interactions, an orthogonality catastrophe analogous to the x-ray edge effect for core levels in a 
metal strongly reduces this singularity to the form (e — A)*^"^", with f3 « 0.3 for typical carbon 
nanotubes. Despite the anomalous tunneling characteristic, good quasiparticles corresponding to 
the unoccupied subband states do exist. 

PACS: 71.10.Pm, 71.20.Tx, 72.80.Rj 

I. INTRODUCTION 

It is well known that one-dimensional (Id) metals are 
proundly influenced by interactions. The generic be- 
havior for a Id system with a single conduction band 
is that of a Luttinger liquid, in which the quasiparticle 
excitations of the non-interacting system are converted 
into spin and charge collective modes. These collective 
modes are largely orthogonal to the bare quasiparticle 
states close to the Fermi surface, hence leading to a dra- 
matic power-law suppression of the tunneling Density of 
States (DOS) at the Fermi energy. A similar vanish- 
ing occurs in the momentum-dependent single-particle 
spectral function, signaling complete breakdown of the 
quasiparticle picture. Indeed, tunneling experiments evi- 
dence the novel nature of the Luttinger liquid much more 
strongly than do other measures such as four-probe re- 
sistance and optical conductivity, which are essentially 
probes of collective modes. There now appears to be in- 
creasing evidence for Luttinger liquid behavior in carbon 
nanotubes, [|lj which are perhaps the most ideal experi- 
mental quantum wires to date. |2,0| 

In this paper, we address the behavior of the tunnel- 
ing DOS far above the Fermi energy. As the voltage 
bias between a tunneling tip and the quantum wire is 
increased, it becomes possible to add an electron not 
only to the conduction band(s), but to higher energy 
unoccupied subbands. Such multiple subband struc- 
ture exists both in carbon nanotubes and semiconduc- 
tor quantum wires. In the former case, the unoccupied 
subbands are essentially states of different angular mo- 
mentum around the graphitic cylinder. For semicon- 
ductor systems, the higher subbands arise from differ- 
ent standing-wave modes transverse to the wire's prop- 
agation direction in the conflnement region. In a non- 
interacting model, the density of states would exhibit van 
Hove singularities at the subband edges. In one dimen- 
sion, these singularities are divergent, giving a contribu- 
tion po{e) "^ ^/m{e — A)~^/-^0(e — A) for energies just 
above the subband edge at e = A (m is the subband 
effective mass). An asymmetric peak structure has in- 
deed been observed in STM measurements of individual 



nanotubes on gold surfaces. 

How do interactions affect these van Hove singulari- 
ties? A simplified, though unphysical model in which the 
mass of the higher subband is taken to be infinite pro- 
vides considerable insight. In this limit the higher energy 
subbands can be replaced by discrete, localized levels. 
The "x-ray edge" problem of a localized level interact- 
ing with a conduction sea was solved by Nozieres and de 
Dominicis ||6|, and is one of the first demonstrations of 
an orthogonality catastrophe. Physically, the core hole 
is "dressed" through interactions with conduction elec- 
trons, which see the hole as a scattering center. A bare 
or undressed hole is then orthogonal to its dressed coun- 
terpart, since an infinite number of conduction electrons 
are available to scatter off of it. This leads to a broaden- 
ing and reduction of the tunneling density of states from 
a sharp delta-function to a power law singularity. While 
this effect is superficially similar to the suppression of 
spectral weight in the Luttinger liquid, it is in fact quite 
distinct. It is due not to the absence of a well-defined core 
hole excitation, but due to its mixing with the conduc- 
tion sea. The distinction is emphasized by the fact that 
x-ray edge singularities are present in any dimension, not 
just in Id. 

Suppose now that the mass of the higher subband is 
taken finite. Since for the tunneling density of states only 
a single particle is being added to the metal, one is faced 
with understanding the behavior of a heavy particle in a 
Luttinger liquid. This problem has been investigated by 
a number of authors in a different context. |Q,|| These 
authors were primarily concerned with the uiobility of 
the heavy particle in response to an external electric field 
at finite temperature. For the tunneling DOS, one is 
interested in a rather different property - essentially, the 
overlap between the two ground states in the presence 
and absence of the heavy particle. This overlap may be 
thought of as a boundary condition changing operator 
g, and is completely outside the Hilbert space of the 
problem in which the heavy particle is always present. 

In this paper, we demonstrate that x-ray edge effects 
persist even for this finite mass case. These reduce the 
naive van Hove singularities in the tunneling DOS by a 



power-law amount. Like the singularities in the original 
x-ray edge problem (but unlike the the Luttinger singu- 
larities at low bias), this modification does not, however, 
signal the destruction of sharp quasiparticle excitations 
in the higher subbands. We argue that a necessary and 
sufficient condition for the presence of such finite-energy 
singularities is a conserved quantum number distinguish- 
ing the states of the higher subband from the conduction 
states. In the case of the carbon nanotube, this is an 
angular momentum quanta. For a semiconductor wire, 
such a good quantum number exists in the ideal case of a 
symmetric confining potential, in which case the second 
subband has odd parity with respect to reflection, while 
the lowest subband has even parity. If such a distinguish- 
ing quantum number is absent, we expect the van Hove 
peak to be rounded and rendered completely nonsingular. 
The case when only forward scattering interactions 
are present between the higher subband particles and 
the conduction sea is asymptotically exactly soluble by 
bosonization methods, as we outline here. As discussed 
in Ref. [|0[ , this forward scattering model is in fact an ex- 
cellent approximation for single- walled carbon nanotubes 
with diameters D > Inm. Within this model, we predict 
a reduced density of states singularity 
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where 0{x) is the heavyside step function, and the or- 
thogonality exponent (3 ~ 0.3 for typical metallic nan- 
otubes. The form of Eq. |l| is expected to apply to the 
second subband in semiconductor quantum wires as well 
(see the concluding remarks for a discussion of experi- 
ments) . 



II. FORWARD SCATTERING CHARGE MODEL 

In this section we present a simple Forward Scatter- 
ing Charge Model (FSCM) describing only forward scat- 
tering interactions, i.e. those processes involving small 
momentum exchange, in the total charge channel. We 
describe the conduction electrons by a Luttinger model, 
valid near the Fermi energy, which we take to be E' = 0, 



i?Q = y~^ / dx vp 
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Here vp is the Fermi velocity, and a =1, | labels the elec- 
tron spin. We have also included an additional "flavor" 
index i — 1 . . . Nf to allow for extra degenerate bands 
(with the same Fermi velocity) at the Fermi energy. For 
metallic carbon nanotubes, Nf = 2 and i should be 
interpreted as a sublattice index. No such special de- 
generacy is present for a semiconductor quantum wire, 
so Nf = 1 in this case. Eq. can be rewritten using 
bosonization. We follow the conventions of Ref. [nO|. One 
has 4'R/L;ia ~ Qi-{'t>ia:±eia) ^ where the dual fields satisfy 



[(j>ia{x),6jp{y)] = -iTT5ij5ap&{x - y) (6(x) is a heavy- 
side step function). Then H^ = J2i a I^^'^oi^ia^^'ia), 
with 
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The slowly varying electronic density in a given channel 
is given by pi^ = ''Pma^Bia+^Lia'pLia = dJia/Tr. Phys- 
ically, 9 can be understood as a displacement or phonon 
field, while cj) carries the phase of the quantum wavefunc- 
tion. 

It is simplest to work in a rotated basis of collec- 



tive modes. For N 
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For Nf = 2, let 0,,p/, = (d^ ± ^^{)|^f2 and 0^± = 
(6*1^ ± 92^)1^/2, with fi = p,(7. With similar definitions 
for the (f> fields, canonical commutators are preserved, 
and Hq —J2alx '^oi^a, (j^a), where a is summed over the 
2A^/ rotated boson fields. 

Because we are interested only in energies near the pu- 
tative van Hove singularity, the unoccupied Id subband 
can be described by a non-relativistic electron operator 
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Here A is the gap to the first subband and m is an effec- 
tive mass. The electron field satisfies {d^{x),dt{x')} — 
daf3S{x — x'). In the case of a carbon nanotube, there 
are actually multiple degenerate subbands at energy A. 
This degeneracy is unimportant within the FSCM, as the 
tunneling DOS involves only states with a single excited 
electron We therefore neglect this additional degeneracy. 
The interactions in the FSCM are written as a single 
term coupling only the total charge density. 



H\„t = 



- dxdx' ptot{x)V{x - x')ptot{x'), (5) 
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Here and in what follows we abbreviate 
the Nf — 2 case. A phenomenological form for the po- 
tential is sufficient for our purposes. We take V{x) = 
exp{—\x\/£^)/y/x^ + W^, modeling the smoothing of the 
interaction on the scale of the wire width by W and in- 
cluding a screening length ^, determined, e.g. by the 
distance to an external gate (any dielectric constant can 
be included by rescaling e^ — *■ e^/e. 

While it is possible to proceed directly with the non- 
local form in Eq. H, near to the van Hove singularity 
(within an energy of order vp/W, up to a weak logarith- 
mic factor) it is sufficient to make the local approxima- 
tion V{x) ^ S{x)Jdx'V{x'). This gives ilmt = L^int, 
with 
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III. SOLUTION OF FSCM 

To determine the effects of the interaction, it is conve- 
nient to employ a path integral formulation. Quantum 
mechanical expectation values are evaluated as functional 
integrals over classical fields in imaginary time with re- 
spect to a measure exp{— JdxdrL), where £ is a La- 
grange density. Standard techniques give 



C — —dx0pdr4>p + d)drd + H, 
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with 



Ti = Hq+Hq+ Hint 
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Here we have defined the plasmon velocity Vp — 
y^vp{vF + {4:Nfe^ /-Kh) \n{^/W)), Luttinger parameter 



("conductance") g = vp/vp, and 7 



= (l-5')- 



It may be tempting to proceed by perturbation the- 
ory in 7. Indeed, for properties of the conduction elec- 
trons at energies small compared to A, this is a perfectly 
sensible procedure: as no heavy particles are present, 
the properties of the conduction sea are completely un- 
affected. However, the same is not true for the tunnel- 
ing density of states. This is obtained from the heavy 
particle Green's function, G{x,t) = {d{x,T)d'' {0,0}), via 
p(e) = T:-^lm Jdk/2TTG{k,iuj -^ uj + iO+). The first per- 
turbative correction to G, obtained from the self-energy 
diagram in Fig. la, is logarithmically divergent. Al- 
though we will not proceed along this route, this log- 
arithmic divergence can be controlled to leading order 
using a renormalization group procedure which treats in 
a self-consistent fashion both this self-energy correction 
and the additional vertex renormalization given by the 
diagram in Fig. lb. The results of this calculation |0] 
are confirmed by a non-perturbative analysis, to which 
we now turn. 




(a) 



(b) 



(c) 



FIG. 1. Diagrams in perturbation theory. A direct attack 
requires a renormalization group resummation of diagrams 
(a) and (b), where the solid line is the heavy-particle propa- 
gator, the dashed line is the Op propagator, and the vertices 
each carry a momentum factor. After the transformation of 
Eqs. thHllI, the residual interactions in Eq. h3 generate only 
irrelevant corrections via diagrams (a) and (c), where now the 
dashed lines represent (pp propagators and the vertices carry 
two momentum factors each. 

The model in Eq. M is solved by a canonical transfor- 
mation, or change of variables in the path integral: 



'p(x) = 9pix)~j / dx'd''ix')dix'), 
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(10) 
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Eqs. |lClfjll| embody the physical process in which the 
conduction sea adiabatically adjusts to the heavy parti- 
cle. In particular, Eq. |l^ represents the depletion of the 
conduction electron density near the heavy particle due 
to Coulomb repulsion. Eq. |l^ represents phase shifts of 
these conduction electrons when the heavy particle is in- 
troduced. Formally, the exponential of the dual (0) field 
in Eq. |l| is a Jordan- Wigner "string" operator which has 
been attached to the heavy particle. 

In the new variables, the hamiltonian density becomes 

n = ^{djpf + ^-^{dx<pf + ni[ld)]+'Hint. (12) 

with a different residual interaction term 
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It might appear that the transformations have actu- 
ally worsened the problem, as the interaction in Eq. 03 
naively appears more complicated than the original in 
Eq. pi However, a closer inspection shows that the new 
couplings in Eq. O arc dimcnsionally weaker by one in- 
verse power of length than the original forms. As the 
original interaction was marginal, the terms in Eq. |l3| 
are in fact irrelevant in the renormalization group sense. 
This can be verified by an explicit calculation of their 
effects on the d Green's function. Apart from a constant 
renormalization of the subband gap, the leading order 
diagrams (Figs. la,lc) give self-energy contributions pro- 
portional to {uj — A)'^. 

The irrelevance of the couplings in Eq. |l^ indicates 
that at long times and distances, the transformed fermion 
and boson correlation functions asymptotically factorize. 
Thus at energies close to the threshold energy A, the d^ 
operator creates good quasiparticles which propagate in- 
dependently of the conduction sea. The tunneling DOS, 
however, involves the addition of a bare electron created 
by the d^ operator. Factorization implies 



G{x,t) = G°(a;,r)/ 
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where G^ is the non-interacting Green's function de- 
scribing free propagation in the unoccupied subband, 
(3 = ^"^ I^Tp'g = (1 — g'^Y /{^Nfg), and A is a momentum 
cutoff (O(fcF)). When Fourier transformed, the space- 
time product above becomes a convolution, which phys- 
ically represents the emission of plasmon waves by the 
added electron. Analytic continuation to real frequency 
gives the modified van Hove singularity in Eq. |l|. 



IV. DISCUSSION 

The preceding analysis demonstrates the persistance 
of a well-defined edge in the tunneling spectrum in the 



presence of forward scattering charge interactions with 
the conduction electrons. The existence of such a finite 
energy singularity hinges on the inability of the heavy 
particle to truly decay or mix with the many other (con- 
duction subband) excitations coexisting at the same en- 
ergy. This is ensured within the FSCM due to heavy 
particle charge conservation. If there are no true distin- 
guishing quantum numbers of the excited state, decay is 
possible and the singularity is rounded, as can be verified 
by explicit calculation using, e.g. Fermi's golden rule. 

We have already argued that for ideal conducting nan- 
otubes and symmetric semiconductor quantum wires, at 
least the first excited subband is protected from decay 
in this way. In any experiment, various non- ideal per- 
turbations will lead to some rounding. Thermal smear- 
ing (from thermal excitation both internally and in the 
tunneling lead) limits the resolution to e > ksT. More 
significant rounding arises from hybridization of the one- 
dimensional electronic states with the bulk. This is prob- 
ably the most significant effect in recent tunneling experi- 
ments using nanotubes on metallic gold substrates. This 
effect could be greatly reduced by using an insulating 
substrate or better, a freely suspended tube. Even for 
nanotubes on an insulating substrate, the asymmetry in 
dielectric constants allows some degree of mixing of dif- 
ferent subband states. Fortunately, this is likely to be 
a weak effect due to the delocalization of the subband 
states around the circumference of the cylinder. Finally, 
impurities or defects lead to smearing of the ideal den- 
sity of states on the scale of the inverse elastic scattering 
time. 

Additional processes left out of the FSCM may lead 
to significant modifications of the physics, although they 
cannot remove the edge singularity. The forward scat- 
tering assumption itself is always correct at energies suf- 
ficiently close to A, since the heavy particle cannot ac- 
comodate a large change in momentum without a corre- 
sponding large increase in its energy. There are, however, 
orward scattering interactions outside the total charge 
channel. The most natural example is the Kondo inter- 
action, Tix = ~^[wr<^^r + '4'l'^'4'l) ' d^'^d. For m = oo, 
one has negligible effects for ferromagnetic (A > 0) cou- 
pling and perfect screening of the heavy spin for antiferro- 
magnetic coupling. For finite m, one can perform a per- 
turbative renormalization group calculation |0] which 
demonstrates that small ferromagnetic coupling remains 
irrelevant for m < oo. Likewise, for antiferromagnetic 
exchange, the problem scales to strong coupling. The na- 
ture of the strong-coupling excitation(s) is not completely 
clear, though preliminary investigations of some lattice 
models suggest the formation of a propagating charge e 
singlet particle. As for the infinite mass model, M we 
expect a large (universal?) contribution to the orthogo- 
nality exponent /3 in this case. In the cases of nanotubes 
and quantum wires, as for most itinerant systems, the 
bare Kondo couplings arc ferromagnetic, and wc thus ex- 
pect insignificant modifications in the spin channel. The 
antiferromagnetic case could perhaps be relevant, how- 
ever, for certain excited states in ladder materials. The 



extra degeneracies in the unoccupied subbands in nan- 
otubes might also lead to some effects of this type for 
pseudospin, but only very close to the band edge as all 
corrections to the FSCM are small in this case. |10| 

We conclude by summarizing the experimental pre- 
dictions for nanotubes and symmetric quantum wires. 
While the tunneling density of states for insulating sys- 
tems should show the bare (e — A)~^/^ van Hove singu- 
larities (up to the two-particle continuum), conducting 
systems should display a reduction of the edge singular- 
ity at the first subband to (e — Ay/^~^. For a lAnm 
diameter nanotube, we estimate /3 ~ 0.3. The singular- 
ities at all higher subbands in metallic quantum wires 
should be rounded. 
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